home *** CD-ROM | disk | FTP | other *** search
/ Loadstar 174 / 174.d81 / star locator (.txt) < prev    next >
Commodore BASIC  |  2022-08-26  |  4KB  |  149 lines
  1. 5 poke55,.:poke56,56:clr
  2. 10 dv=peek(186):ifdv<8thendv=8
  3. 15 poke53280,0:poke53281,0:print"[147]"
  4. 16 poke53371,0
  5. 25 poke53272,31
  6. 30 ad=49152
  7. 35 sysad:sysad+12
  8. 100 rem rigorous precession
  9. 105 rem (j2000.0)
  10. 110 p1=(NULL):r1=p1/180
  11. 112 bx$="[158]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_"
  12. 113 sb$="[153]_[157]_[157]_[157]_[157]_[157]_[157]_"
  13. 114 tq=ti
  14. 115 e0=2000:a$="a"
  15. 120 print"[147]":sysad+9,6:poke53272,31
  16. 121 print"[153]   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
  17. 122 printtab(3)sb$;""tab(33)sb$
  18. 123 print"[153]   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
  19. 124 poke214,2:print:printtab(9)"[129][211][212][193][210] [204][207][195][193][212][207][210] [208][210][207][199][210][193][205]!"
  20. 125 poke214,4:print:printtab(5)"[158][201]nitial [197]poch [155](yr) ";:l9%=4:gosub670:ie=q9
  21. 130 print:printtab(5)"[158][198]inal [197]poch [155](yr):";:l9%=7:gosub670:fe=q9
  22. 135 ifabs(ie-fe)<=500 then150
  23. 140 poke214,10:print:printtab(9)"[212]ime span too long!"
  24. 141 tq=ti
  25. 142 ifti<tq+200then142
  26. 145 goto120
  27. 150 print"[147]":sysad+9,17:gosub415
  28. 151 print"[158]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
  29. 152 printbx$;""tab(38)bx$
  30. 153 print"[158]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
  31. 155 ifa$="n"then185
  32. 160 print"":printtab(9)"[153] [201]nitial [197]poch [146][158]";ie
  33. 165 printtab(3)"[154][210].[193]. ";
  34. 170 gosub350:a=v
  35. 175 printtab(3)"[150][196]ec  ";
  36. 180 gosub350:d=v
  37. 185 ma=0:md=0:gosub570
  38. 190 a2=a1:d2=d1
  39. 195 ifa$="n"then220
  40. 200 printtab(9)"[153] [208]roper [205]otion in "
  41. 205 printtab(3)"[159][210].[193]. [155](sec/yr) ";:l9%=7:gosub670:m1=q9
  42. 210 printtab(3)"[159][196]ec. [155](arcsec/yr) ";:l9%=7:gosub670:m2=q9
  43. 215 printtab(1)"[158]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^":sysad+9,8
  44. 220 ma=m1:md=m2:gosub570
  45. 225 v=a1:gosub380
  46. 230 printtab(9)"[153] [198]inal [197]poch [146][158] ";fe
  47. 235 printtab(3)"[159][210].[193]. [155](h,m,s):[158] ";
  48. 240 printv1;v2;v3
  49. 245 v=d1:gosub380
  50. 250 printtab(3)"[159][196]ec. [155](d,m,s):[158] ";s$;
  51. 255 printv1;v2;v3
  52. 260 a3=(a1-a2)*3600/ny
  53. 265 d3=(d1-d2)*3600/ny
  54. 270 printtab(9)"[153] [208]roper [205]otion in "
  55. 275 printtab(3)"[159][210].[193]. [155](sec/yr):[158] ";
  56. 280 printa3
  57. 285 printtab(3)"[159][196]ec. [155](arcsec/yr):[158] ";
  58. 290 printd3
  59. 295 print
  60. 300 printtab(9)"[159] [211]elect [207]ne [208]lease:[146] "
  61. 305 printtab(3)"[150]([193])[153]nother [211]tar   [150]([206])[153]ew [198]inal [197]poch"
  62. 315 printtab(11)"[150]([213])[153]niverse [205]enu
  63. 320 [161]a$:[139]a$[179][177]"a"[175]a$[179][177]"n"[175]a$[179][177]"u"[167]320
  64. 325 [139]a$[178]"a"[167]150
  65. 330 [139]a$[178]"n"[167]110
  66. 335 [139]a$[178]"u"[167]3050
  67. 350 [143] input sexagesimal
  68. 355 s[178]1:[153]"conthour/defdegree (def+/-):";:l9%[178]3:[141]670:q9[178][197](q9$):v$[178]q9$
  69. 356 [153][163]9)"contminutes:";:l9%[178]4:[141]670:v2[178]q9
  70. 357 [153][163]9)"contseconds:";:l9%[178]3:[141]670:v3[178]q9
  71. 360 [139][200](v$,1)[178]"-"[167]s[178][171]1
  72. 365 v1[178][182]([197](v$))
  73. 370 v[178]s[172](v1[170]v2[173]60[170]v3[173]3600)
  74. 375 [142]
  75. 380 [143] output sexagesimal
  76. 385 s$[178]"+":[139]v[179]0[167]s$[178]"-"
  77. 390 v[178][182](v):v1[178][181](v)
  78. 395 vm[178]60[172](v[171]v1):v2[178][181](vm)
  79. 400 v3[178]60[172](vm[171]v2)
  80. 405 [142]
  81. 410 :
  82. 415 [143] pression parameters
  83. 420 ny[178]fe[171]ie
  84. 425 t0[178](ie[171]e0)[173]100:t1[178]ny[173]100
  85. 430 t2[178]t1[172]t1:t3[178]t1[172]t1[172]t1
  86. 435 h1[178]2306.2181:h2[178]1.39656
  87. 440 h3[178].000139:h4[178].30188
  88. 445 h5[178].000345:h6[178].017998
  89. 450 k1[178]1.09468:k2[178].000066
  90. 455 k3[178].018203
  91. 460 l1[178]2004.3109:l2[178][171].8533
  92. 465 l3[178][171].000217:l4[178][171].42665
  93. 470 l5[178][171].000217:l6[178][171].041833
  94. 475 w[178](h1[170]h2[172]t0[170]h3[172]t0[172]t0)[172]t1
  95. 480 zt[178]w[170](h4[170]h5[172]t0)[172]t2[170]h6[172]t3
  96. 485 zd[178]w[170](k1[170]k2[172]t0)[172]t2[170]k3[172]t3
  97. 490 th[178](l1[170]l2[172]t0[170]l3[172]t0[172]t0)[172]t1
  98. 495 th[178]th[170](l4[170]l5[172]t0)[172]t2[170]l6[172]t3
  99. 500 zt[178]zt[172]r1[173]3600:zd[178]zd[172]r1[173]3600
  100. 505 th[178]th[172]r1[173]3600
  101. 510 [143] zt,zd,th = euler amgles
  102. 515 :
  103. 520 [143] rotation matrix
  104. 525 s1[178][191](zt):c1[178][190](zt)
  105. 530 s2[178][191](zd):c2[178][190](zd)
  106. 535 s3[178][191](th):c3[178][190](th)
  107. 540 xx[178]c1[172]c3[172]c2[171]s1[172]s2
  108. 545 yx[178][171]s1[172]c3[172]c2[171]c1[172]s2:zx[178][171]s3[172]c2
  109. 550 xy[178]c1[172]c3[172]s2[170]s1[172]c2
  110. 555 yy[178][171]s1[172]c3[172]s2[170]c1[172]c2:zy[178][171]s3[172]s2
  111. 560 xz[178]c1[172]s3:yz[178]s1[172]s3:zz[178]c3
  112. 565 [142]
  113. 570 [143] proper motion correction
  114. 575 a0[178](a[170]ma[172]ny[173]3600)[172]15[172]r1
  115. 580 d0[178](d[170]md[172]ny[173]3600)[172]r1
  116. 585 :
  117. 590 [143] sherical rectangular
  118. 595 sa[178][191](a0):ca[178][190](a0)
  119. 600 sd[178][191](d0):cd[178][190](d0)
  120. 605 x0[178]ca[172]cd:y0[178]sa[172]cd:z0[178]sd
  121. 610 [143] 3-d transformation
  122. 615 x1[178]x0[172]xx[170]y0[172]yx[170]z0[172]zx
  123. 620 y1[178]x0[172]xy[170]y0[172]yy[170]z0[172]zy
  124. 625 z1[178]x0[172]xz[170]y0[172]yz[170]z0[172]zz
  125. 630 [143] rectangular sperical
  126. 635 a1[178][193](y1[173]x1)
  127. 640 [139]x1[179]0[167]a1[178]a1[170]p1
  128. 645 [139]a1[179]0[167]a1[178]a1[170]2[172]p1
  129. 650 a1[178]a1[173](r1[172]15):[143] final r.a.
  130. 655 d1[178][193](z1[173][186](x1[172]x1[170]y1[172]y1))
  131. 660 d1[178]d1[173]r1:[143] final dec.
  132. 665 [142]
  133. 670 q9$[178]"":[151]198,.
  134. 675 [161]a$
  135. 680 [151]646,[187](1)[172]15[170]1:[153]"*cmd";:[139]a$[178]""[167]675
  136. 685 [139]a$[178][199](13)[167][153]" ":q9[178][197](q9$):[142]
  137. 690 [139](a$[178][199](20)[175][195](q9$))[167]q9$[178][200](q9$,[195](q9$)[171]1):[137]720
  138. 695 [139][195](q9$)[177][178]l9%[167][158]ad[170]9,1:[137]675
  139. 700 [139](a$[177][178]"0"[175]a$[179][178]"9")[176]a$[178]"."[176]a$[178]"-"[176]a$[178]"+"[167]710
  140. 705 [137]675
  141. 710 q9$[178]q9$[170]a$
  142. 715 [153]"sys"a$;:[158]ad[170]9,7:[137]675
  143. 720 [153]" cmdcmd cmd";:[137]675
  144. 3050 [158]ad[170]15
  145. 3060 [153]"loadstopload"[199](34)"b.universe iii"[199](34)","dv
  146. 3070 [153]"run28"
  147. 3080 [151]631,13:[151]632,13:[151]198,2:[128]
  148. 10000 d[178][194](186):n$[178]"0:star locator":[159]15,d,15,"s"[170]n$:[160]15:[148]n$,d:[128]
  149.